Strongly Interacting Polaritons in Coupled Arrays of Cavities 
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The experimental observation of quantum phenomena in strongly correlated many particle systems 
is difficult because of the short length- and timescales involved. Obtaining at the same time detailed 
control of individual constituents appears even more challenging and thus to date inhibits employing 
such systems as quantum computing devices. Substantial progress to overcome these problems has 
been achieved with cold atoms in optical lattices, where a detailed control of collective properties 
is feasible but it is very difficult to address and hence control or measure individual sites. Here 
we show, that polaritons, combined atom and photon excitations, in an array of cavities such as 
a photonic crystal or coupled toroidal micro-cavities, can form a strongly interacting many body 
system, where individual particles can be controlled and measured. All individual building blocks 
of the proposed setting have already been experimentally realised, thus demonstrating the potential 
of this device as a quantum simulator. With the possibility to create attractive on-site potentials 
the scheme allows for the creation of highly entangled states and a phase with particles much more 
delocalised than in superfluids. 

PACS numbers: 03.67.Mn, 42.50.Dv, 73.43.Nq, 03.67.-a 



INTRODUCTION 

Physical systems where the interactions between the 
constituent particles are weak and can be treated per- 
turbatively are mostly well understood. A precise un- 
derstanding of strongly correlated systems however ap- 
pears to be much harder to obtain because the experi- 
mental observation of several important mechanisms is 
less feasible. Many processes take place at high energy 
and therefore on fast timescales while precise control over 
individual constituent particles is very hard to achieve. 
This unattainability of detailed control lies at the heart 
of the difficulties to experimentally implement quantum 
information processes. 

Nonetheless, substantial progress has been achieved re- 
cently by considering cold atoms trapped in an optical 
lattice which can effectively be described by a Bose- 
Hubbard Hamiltonian. This idea lead to several impor- 
tant experiments emulating condensed matter phenom- 
ena. The most prominent example is probably the obser- 
vation of the superfluid to Mott insulator quantum phase 
transition Q in a three dimensional lattice, but also other 
geometries like a Kagome lattice Q or a Tonks Girardeau 
gas Q are being studied. Even theoretical ideas for the 
implementation of quantum information processing have 
been pursued by entangling neighbouring atoms via con- 
trolled collisions 0, Q . 

In addition, these systems could act as quantum simu- 
lators, where the dynamics of a rather complex Hamilto- 
nian can be simulated in an experiment by eather creat- 
ing the Hamiltonian in question, or Trotter decomposing 
the evolution into small steps, where the dynamics of 
each step is generated by a much simpler Hamiltonian 



0i. 

Despite their success, optical lattices are in all their 
applications limited by an intrinsic draw back: It is ex- 
ceedingly difficult to access individual lattice sites, since 
their separation is only half of the employed optical wave- 
length. 

Here we propose an alternative way to create an effec- 
tive Bose-Hubbard model and other Hamiltonians, which 
does not suffer from this problem. We consider two re- 
alisations of an array of cavities and study the dynamics 
of polaritons, combined atom photon excitations, in this 
arrangement. Since the distance between adjacent cavi- 
ties is considerably larger than the optical wavelength of 
the resonant mode, individual cavities can be addressed. 

Photon hopping occurs between neighboring cavities 
while the repulsive force between two polaritons occupy- 
ing the same site is generated by a large Kerr nonlinearity 
that occurs if atoms with a specific level structure usually 
considered in Electromagnetically Induced Transparency 



Id, lUl, yj, liJ, yj, lis! interact with light. Each cav- 



ity is interacting with an ensemble of these atoms, which 
are driven by an external laser, see fig. [TJ By varying 
the intensity of the driving laser, the nonlinearity can 
be tuned and hence the system can be driven through 
the superfluid to Mott insulator transition. In partic- 
ular, the driving laser may be adjusted for each cavity 
individually allowing for a much wider range of tuning 
possibilities than in an optical lattice, including attrac- 
tive interactions. 

Our system can be described in terms of three species 
of polaritons, where under the conditions specified below, 
the species which is least vulnerable to decay processes 
is governed by the effective Bose-Hubbard (BH) Hamil- 
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FIG. 1: Our model consists of an array of cavities, where 
photon hopping occurs due to the overlap (shaded green) of 
the light modes (green lines) of adjacent cavities. Atoms in 
each cavity (brown) , which are driven by external lasers (blue) 
give rise to an on site repulsion. 



tonian, 

R {R,R') 

creates a polariton in the cavity at site R and the 
parameters k and J describe on site repulsion and inter 
cavity hopping respectively. 

The most promising candidates for an experimental 
realisation are toroidal or spherical micro-cavities, which 
are coupled via tapered optical fibres These cavi- 

ties can be produced and positioned with high precision 
and in large numbers. They have a very large Q-factor 
(> 10^) for light that is trapped as whispering gallery 
modes and efficient coupling to optical fibres [17| as well 
as coupling to Cs-atoms in close proximity to the cavity 
via the evanescent field l^, 1^ have been demonstrated 
experimentally very recently. 

In the longer term, photonic crystals represent an ap- 
pealing alternative as they offer the possibility for the 
fabrication _of large arrays of cavities in lattices or net- 
works |2C 
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|22| . |23|. Obviously, our concept is not 
limited to certain geometries, as it is the case for the op- 
tical lattice. Any arrangement of the cavity array may 
be considered. 

We first derive the Hamiltonian ([T|) for the considered 
structure and then present a theoretical analysis of the 
feasible parameter range, which is backed up by a full nu- 
merical demonstration of the supcrfluid to Mott insulator 
transition including experimental imperfections. 



sorption processes, and limit our considerations to lin- 
ear, isotropic dielectric media. The electromagnetic field 
may be represented by a vector potential A and a scalar 
potential <!> which obey the gauge conditions, $ = and 
V • {e{f) A) ~ [23|. A can be expanded in Wannier 
functions, w^, each localised at one single cavity at lo- 
cation R. We describe this single cavity by the dielectric 
function ej^{r) such that the Wannier functions satisfy 
the eigenvalue equation. 



V X (V X uJj^j 



0, 



(3) 



where the eigenvalue uJq is the square the resonance fre- 
quency of the cavity which is independent of R due to 
the periodicity. We assume that the Wannier functions 
decay strongly enough outside the cavity such that only 
Wannier modes of nearest neighbour cavities have non- 
vanishing overlap. 

In terms of the creation and annihilation operators of 
the Wannier modes, a^-, and a^, the Hamiltonian of the 
field can be written, 



E 



h.c. 



{R,R') 



(4) 

Here sum of all pairs of cavities which 

are nearest neighbours of each other. Since a <^ 1, we 
neglected rotating terms which contain products of two 
creation or two annihilation o pera tors of Wannier modes 
in deriving ([4]). a is given by 23, 2^, 



(fr (e^(f) - e(f)) w^w^, 



\R-R'\ ^R, (5) 



and can be obtained numerically for specific models [26| . 

The model (|4]) provides an excellent approximation to 
many relevant implementations such as coupled photonic 
crystal micro-cavities or fibre coupled toroidal micro- 
cavities. Furthermore, it allows for the observation of 
state transfer (27| and entanglement dynamics and prop- 
agation for Gaussian states |28l . [29j as the Hamiltonian 
is harmonic. In the next section we present a possible 
realisation of the repulsive term in the Hamiltonian ([1]). 



THE EFFECTIVE ON SITE REPULSION 



QUANTISED ELECTROMAGNETIC FIELD IN A 
PERIODIC ARRAY OF CAVITIES 

We consider a periodic array of cavities, which we de- 
scribe here by a periodic dielectric constant, 

e(f) = e(f+i?); R^nR, (2) 

for all tupels of integers ft = {nx,ny^nz). We assume 
the dielectric constant e to be real, i.e. we neglect ab- 



To generate a repulsion between polaritons that are lo- 
cated in the same cavity, we fill the cavity with 4 level 
atoms of a particular level structure that are driven with 
an external laser in the same manner as in Electromag- 
netically Induced Transparency, see figure [S] The transi- 
tions between levels 2 and 3 are coupled to the laser field 
and the transitions between levels 2-4 and 1-3 couple via 
dipole moments to the cavity resonance mode. It has 
been shown by Imamoglu and co-workers, that this atom 



3 



A- 




5i3 yijjc 



of atoms is large, A'' ^ 1, the Hamiltonian ([6]) truncated 
to the subspace of at most two excitations can be diag- 
onahsed. Let us therefore define the following creation 
(and annihilation) operators: 



pI = 



A{A + 5) 



A{A - 5) 



A 



ga' 



2 

A-5 



-S- 



13 



^13 



(7) 



FIG. 2: The level structure and the possible transitions of 
one atom, cue is the frequency of the cavity mode, SIl is the 
Rabi frequency of the driving by the laser, gis and g24 are the 
parameters of the respective dipole couplings and S, A and e 
are detunings. 



cavity system can show a very large nonlinearity 11| , and 
a similar nonlinearity has recently been observed experi- 
mentally [3^ . 

In a rotating frame with respect to Hq = 

LUC (a'^a + ^) + Ei^i (^00-22 + + ^LOccri^^ 

Hamiltonian of the atoms in the cavity reads, 



the 



N 

E 

N 

E 

i=i 



(A + e)ai 



913 CT^3 a'f 



924 CT^4 



h.c 



(6) 



where a^i = \kj){lj \ projects level I of atom j to level k of 
the same atom, cue is the frequency of the cavity mode, 
r^L is the Rabi frequency of the driving by the laser and 
(713 and (724 are the parameters of the dipole coupling 
of the cavity mode to the respective atomic transitions 
which are all assumed to be real. 

Since we are only interested in the situation, where 
there is on average one excitation in each cavity and since 
the probability to have two excitations in one cavity is 
suppressed by the Kerr nonlinearity, it is appropriate to 
only consider states with at most two excitations per cav- 
ity. All atoms interact in the same way with the cavity 
mode and hence the only relevant states are Dicke type 
dressed states. Neglecting two photon detuning and cou- 
pling to level 4, the Hamiltonian ^ can be written as a 
model of three species of polaritons in this subspace. 



Polaritons 

In the case where (724 = and e = , level 4 of the 
atoms decouples from the rest of the two excitation man- 
ifold ij, [3IJ . If we furthermore assume that the number 



where g = VNgi3, B = ^/g^ 
^ii and 



•^12 



A ^ V4.B^ + S^, 



Eiv 7 
, = 1^31- 



The operators pj, and p[_ describe polaritons. quasi 
particles formed by combinations of atom and photon 
excitations. By looking at their matrix representation in 
the subspace of at most two excitations, one can see that 
in the limit of large atom numbers, ^ 1, they satisfy 
bosonic commutation relations, 



IP J, Pi. 



= 



Pj^Pi 



J, 1 = 0,+,-. (8) 



Pq, and pi thus describe independent bosonic parti- 
cles. In terms of these polaritons, the Hamiltonian ^ 
for 324 = and £ = reads, 



[Hi]g^^^o,e=o MoPoPo + ^^+P+P+ + A*- 



P-P- 



(9) 



where the frequencies are given by /iq = 0, /i+ = {S — 
A)/2 and ^_ = (J + A)/2. We will now focus on the 
dark state polaritons, pj, which can be described by the 
Hamiltonian ([T]) as we shall see. 



Perturbations 



To write the full Hamiltonian Hj, in the po- 
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and 



lariton basis, we express the operators 

J2jLi '^24 ™ terms of p^, p\^ and p^_. 
In the subspace of at most two excitations, the coupling 
of the polaritons to the level 4 of the atoms via the dipole 
moment 1724 reads. 




-.924 



9^1 



Sl^p'o +h.c. 



where Sl^^ = -^LjLif^ii- In deriving 



(10) 
we made 

use of the rotating wave approximation: In a frame rotat- 
ing with respect to the polaritonic creation operators 
rotate with the frequencies /iq, fJ-+ and fi-. Furthermore, 
the operator 5*^4 rotates at the frequency 2fio, (c.f. 
Hence, provided that 



\g24\, |e|, |A| < 



(11) 
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all terms that rotate at frequencies — {n+ ± or 
jio ± M+ or ± can be neglected, which eliminates 
all interactions that would couple and to the re- 
maining polariton species. 

For (724 ^ |A|, the coupling to level 4 can be treated 
in a perturbative way. This results in an energy shift of 
2 K with 



Ng 



13 ' 



(12) 



and in an occupation probability of the state of one 
excitation of -'2k/ A, which will determine an effective 
decay rate for the polariton via spontaneous emission 
from level 4 (see below). Note that k > for A < and 
vice versa. In a similar way, the two photon detuning 
e leads to and energy shift of e for the polariton 

Pq, which plays the role of a chemical potential in the 
effective Hamiltonian. 

Hence, provided pip holds, the Hamiltonian for the 
dark state polariton p\ can be written as 



2 2 



(13) 



in the rotating frame. Next we turn to estimate the ef- 
fects of loss and decoherence. 



THE COMPLETE PICTURE 

Now we look at an array of cavities which forms a 
periodic dielectric medium, see fig. [TJ The first term of 
the Hamiltonian ^ has already been incorporated in the 
polariton analysis for one individual cavity. The second 
term transforms into the polariton picture via ([7]). To 
distinguish between the dark state polaritons in differ- 
ent cavities, we introduce the notation to label the 

polariton Pq in the cavity at position R. We get. 



(15) 



"terms for other polariton species" 



where we have applied a rotating wave approximation. 
Contributions of different polaritons decouple due to the 
separation of their frequencies /io, ^J■+ and As a con- 
sequence the Hamiltonian for the polaritons p^, takes on 
the form ([1]), with 



J = 



2^c»i 



a, 



(16) 



where we have assumed a negligible two photon detuning, 
e« 0. 



Spontaneous emission and cavity decay 



NUMERICAL ANALYSIS 



Since level 2 of the atoms is metastable and hence its 
decay rate negligible on the relevant time scales, the dark 
state polaritons Pq suffer loss and decoherence via two 
processes only: The photons are subject to cavity decay 
at a rate Tc and spontaneous emission from the atomic 
level 4 at a rate r4 induces an effective decay via its small 
but nonvanishing occupation. 

Taking into account the definition ([7| of p^ and the 
effective decay of Pq induced by the occupation of level 
4, we obtain the following decay rate T for the polaritons 



ph 



r = 



2nl 



Ng 



c 



g^Ngh 



{Ngj, 



(14) 



Due to (724 <C |A|, spontaneous emission from level 4 
affects the polaritons much less than cavity decay. 

For successfully observing the dynamics and phases of 
the effective Hamiltonian (fT3|) , the repulsion term n needs 
to be much larger than the damping F. Choosing J^l <C 
^fN 313, one can achieve a ratio of k/F ^ 0.4 for photonic 
band gap cavities (s^ . while for toroidal micro-cavities 
the much larger value of k/T ^ 1.1 x lO'^ is possible, 
which makes them an ideal candidate for an experimental 
implementation 17, 18, S^]. In the next section we turn 
to combine the on site repulsion for polaritons with the 
photon hopping between neighbouring cavities. 



To provide evidence for the accuracy of our approach, 
we present a numerical simulation of the dynamics of one 
dark state polariton in an array of three toroidal micro- 
cavities with periodic boundary conditions. Parameters 
from [33, 513 = 2.5 x lO^s-i, (?24 = 2.5 x lO^s"!, F4 = 
1.6 X 10'^s-\ Tc = 0.4 X 10^s~\ and assuming 5=l.Qx 
lO^^s-i, A = -1.0 X 10"s-\ N = 10000, e = 0, 2wca = 
0.4 X 10*^s"i and VIl = 2.5 x 10"s~i yields values of 
K = 1.6 X lOV^ and J 2.0 x lOV^ 

The polariton is initially located in cavity 2, and we 
compare its dynamics to that of a pure BH model. In 
doing so we look at the number of polaritons in one cav- 
ity/site, ni ~ {p\pi) for site Z, and the number fluctua- 
tions, A; = {{pjpi)'^) — {pjpi)'^. Figure [3] shows cavity 
polariton numbers and number fluctuations for cavities 
1 and 2, in the upper plot while the lower plot shows the 
differences between the cavity polariton number in one 
cavity and the excitation number at one site for the pure 
BH model, Sni = [ri/] cavities ~ ["'Ibh' ^^^^ ^he dif- 
ference between both models in the number fluctuations. 



SAi = [A,] 



cavities 



-[A/] 



BH- 



All our simulation results 
show a single trajectory of a quantum jump simulation 
[s^l. Within the simulated time range the probability of 
a decay event remained less than 2 x 10^'^. 

A key question is of course, whether the polaritons 
undergo the Mott insulator to supcrfluid phase transition 



FIG. 3: Upper plot: Cavity polariton numbers and number 
fluctuations for cavities 1 and 2. Lower plot: Difi^erences be- 
tween the cavity polariton number in one cavity and the ex- 
citation number at one site for the pure BH model, Srii, and 
difference between both models in the number fluctuations, 
SAi for cavities/sites 1 and 2. Deviations from the pure BH 
model are less than 10""^. 

in the same way as a pure BH model. We give evidence 
for that in the next section. 



The phase transition 

An important feature of the Bose-Hubbard Hamiho- 
nian ([T|) arc two phases of the matter it describes: One 
is a supcrfluid phase, where the hopping term dominates 
over the on site repulsion, J 2> k, and the polaritons 
pj){R) move almost freely between the cavities. As a con- 
sequence one can observe large fluctuations in the number 
of polaritons in each individual cavity. The other phase is 
the Mott insulator, in which the on site repulsion is much 
stronger than the hopping term J <^ k and double occu- 
pancies of a cavity are strongly suppressed. Provided the 
system contains on average one polariton per cavity, this 
implies that in the Mott phase the number of polaritons 
is exactly one for each cavity and number fluctuations 
are very small. 

A numerical simulation of this phase transition is 
shown in figure 3) Again we consider the parameters 
for coupled toroidal micro-cavities, but now we assume 
that the driving laser is tuned from very weak driving to 
strong driving throughout the experiment. We look at 




0.2 0.4 0.6 0.8 1.0 

t in 10-6 g 



FIG. 4: Log-plot of k (blue) and J (green) as a function of 
time, upper plot, number fluctuations of cavity polaritons in 
cavity 1, Ai, as a function of time, middle plot and difference 
between number fluctuations in cavity polaritons and number 
fluctuations in a pure BH model, SAi, as a function of time, 
lower plot. Parameters are the same as in figure [3] except for 
N = 1000, A = -0.2 X 10"s"^ and 2ujca = 0.12 x 10**s"\ 
Deviations from the pure BH model are about 5%. 

four coupled cavities with periodic boundary conditions, 
which are initially prepared such that there is one polari- 
ton Pq in each. Figure |3] shows how the values of k and J 
vary in time due to the tuning (upper plot), the number 
fluctuations of cavity polaritons in cavity 1, Ai (middle 
plot), and the difference between number fluctuations in 
cavity polaritons and number fluctuations in a pure BH 
model, 5Ai (lower plot). The expectation value for the 
number of polaritons, {pqPo) would be equal to one for 
an exact BH model. This also holds for the two cavities 
with deviations less than 3 x 10^^. The probability for 
a decay event is about 15% for the total simulated time 
range. Next we mention some applications of the setup 
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0-2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

t in 10^^ s 

FIG. 5: Overlap with the momentary ground state Ogg = 
|((/)(t)|0ga(t))|, overlap with the Wn state ow = |(0(t)|W'jv)| 
and the number fluctuations Ai for the efi^ective Hamiltonian 
HI) with 4 cavities where J = lO^s"^ and k is linearly tuned 

form Kmin = — 2 X lO^S^^ to Kmax = — 4 X lO^S^""^. 

when operated with attractive on-site potential. 

Attractive on-site potential 

As mentioned above, if the detuning A is positive, the 
on-site potential becomes attractive. The gromid state of 
the system with polaritons on N sites and a strong at- 
tractive on-site potential is then a highly entangled Wn 
state; \Wn) = :^SjLil-^i)j where \Nj) is the state 
with all TV polaritons on site j [s^l • Although the ground 
state becomes degenerate in the limit k — oo, one can 
generate the Wn state with high fidelity since the hop- 
ping terms do not induce transitions between the nearly 
degenerate low energy states. The procedure is robust 
against disorder as long as the hopping J is stronger than 
the fluctuations of k. In the Wn state regime, the polari- 
tons enter a physically different phase, where the number 
fluctuations at one site become A, = A^ — 1. 

Figure [5] shows results of a simulation for the effective 
Hamiltonian ([T|) with 4 cavities where J = 10^s~^ and 
K is tuned form Kmin = — 2 x 10^s~^ to K^ax = — 4 x 
10 ''s"-^. The plotted quantities are the overlap with the 
momentary ground state |(<^(t)|(/)gs(t))|, the overlap with 
the Wn state |((/)(t)|VIOv)| and the number fluctuations 
Ai. 

Switching between positive and negative on-site poten- 
tial is possible since the detuning A can be varied via the 
Stark shift induced by an applied electric field. We now 
turn to discuss how to distinguish different phases in a 
measurement. 



MEASURING THE STATE OF THE 
POLARITONS 

A possible way to detect whether or not there is a 
polariton present in a certain cavity is to make use of 



resonance fluorescence. To that end an additional probe 
laser is resonantely coupled to the transition between an 
atomic level 5, which is well separated from the levels 
1-4, and the level 2. The intensity of the probe laser 
needs to be weak enough such that its Rabi frequency f2p 
fulfils ilp ^ — /io|, l/i- — /io|). In that case the probe 
laser couples a polariton to the extra level 5 while 
its coupling to the other polariton species is negligible. 
Hence whenever resonant light is detected, at least one 
polariton pj, must have been present in the cavity. 

Provided there is on average one polariton per cavity, 
repeated measurements at several cavities reveal whether 
there has always been a polariton in each cavity (Mott 
phase) or whether there are every now and then empty 
cavities indicating polariton number fluctuations (super- 
fluid phase). 

Furthermore the fluorescence spectrum contains infor- 
mation about the energy differences between level 5 and 
the eigenstates of the polaritonic BH model. In the limit 
of negligible Rabi frequency of the probe laser and negli- 
gible linewidth of level 5, the exact spectrum of ([T]), which 
is gapped in the Mott phase and gapless in the superfiuid 
phase, could be observed. In the real situation, the finite 
Rabi frequency of the probe laser has to be taken into 
account and the linewidth of level 5 limits the accuracy 
of the measurement. 

CONCLUSIONS 

We have shown, that polaritons, combined atom - pho- 
ton excitations, can form an effective quantum many 
particle system, which is described by a Bose-Hubbard 
Hamiltonian, and numerically demonstrated that an ob- 
servation of the Mott insulator to superfiuid phase tran- 
sition is indeed feasible for experimentally realisable pa- 
rameters. 

In contrast to earlier realisations in optical lattices, our 
scenario has the advantage that single sites of the BH 
model can be addressed individually. This feature opens 
up various new experimental possibilities: One can for 
example operate one part of the system in the Mott insu- 
lator regime while the rest is in the superfiuid phase and 
study the physics that happens at the boundary. This 
boundary may even be moved around during the exper- 
iment. Local properties of the system may be probed 
[35I . [sl] . Moreover being able to address individual sites 
is a prerequisite for the implementation of quantum in- 
formation processes. Hence we expect our propos al to 
lead to significant progress in that direction [371, l38| . 

In addition, our setup offers the possibility to create a 
BH type Hamiltonian with an attractive on site poten- 
tial, which is readily achieved by choosing the detuning 
A to be positive, and may be used to generate highly en- 
tangled states [39|. Even designing many body systems 
with long range interactions seems feasible since toroidal 
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micro-cavities are usually coupled to light in optical fi- 
bres This fibre can the couple many cavities in a 
row without a significant decrease of the light field along 
it. 
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